A scalar field dark matter model and its role in the large scale structure formation in 
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In this work we present a model of dark matter based on scalar-tensor theory of gravity. With 
this scalar field dark matter model we study the non-linear evolution of the large scale structures 
in the universe. The equations that govern the evolution of the scale factor of the universe are 
derived together with the appropriate Newtonian equations to follow the non-linear evolution of the 
structures. Results are given in terms of the power spectrum that gives quantitative information on 
the large-scale structure formation. The initial conditions we have used are consistent with the so 
called concordance ACDM model. 



I. INTRODUCTION 

The standard model of cosmology is supported by three 
main astronomical observations: the surveys of super- 
novae la, the cosmic microwave background radiation 
(CMB), and the primordial nucleosynthesis. These ob- 
servations together with other modern cosmological ob- 
servations, like galaxies surveys (SDSS, 2dF), galaxy ro- 
tation curves, the Bullet Cluster observation, studies of 
clusters of galaxies, establish that the universe behaves as 
dominated by dark matter (DM) and dark energy. How- 
ever, the direct evidence for the existence of these invis- 
ible components remains lacking. Several theories that 
would modify our understanding of gravity have been 
proposed in order to explain the large scale structure for- 
mation in the universe and the galactic dynamics. The 
best model we have to explain the observations is the 
ACDM model, i.e., the model of cold dark matter (CDM) 
-non-relativistic particles of unknown origin- with cos- 
mological constant (A), in particular, this model explains 
very well the universe on scales of galaxy clusters and up 



The ACDM model has become the theoretical 
paradigm leading the models of the universe to explain 
the large scale structure (LSS) formation and several 
other observations. Where "large" means scales larger 
than 1 Mpc -about the size of the group of galaxies that 
our Milky Way belongs. Together with the cosmic infla- 
tion theory, this model makes a clear prediction about the 
necessary initial conditions that the universe has to have 
in order to have the structures we observe and that those 
structures build hierarchically due to a gravitational in- 
stability. One of its main predictions is that the density 
profile of galaxies, clusters of galaxies, and so on, is of 
the form [2,3, 



PNFw{r) = 



Po 



(r/ro)(l4-r/ro)2 



a density profile known as Navarro-Frenk- White profile 
(NFW). Parameters po and tq must be fitted, for exam- 
ple, using rotation curves of galaxies. 

The ACDM model and its success in explaining several 
observations -this is why this model is also known as the 



concordance model- yields the following conclusions: On 
large scales, the universe is homogeneous and isotropic, as 
described by the Friedmann-Lamaitre- Robertson- Walker 
(FLRW) metric. The geometry of the universe is flat, 
as predicted by inflation. The dark matter is cold (non- 
relativistic at decoupling epoch) . The initial density fluc- 
tuations were small and described by a Gaussian random 
field. The initial power spectrum of the density fluctu- 
ations was approximately the Harrison-Zeldovich spec- 
trum (P(fc) cx fc", n = 1) d [i]. 

In terms of the composition of the universe, the above 
conclusions can be summarized as follows: Bubble's con- 
stant (Expansion rate of the universe at the present 
epoch): Hq = 73.2 ± 3.1 km/s/Mpc. Density param- 
eter (combined mass density of all kind of mass and 
energy in the universe, divided by the critical density): 
ilo = 1.02 ± 0.02. Matter density parameter (combined 
mass density of all forms of matter in the universe, di- 
vided by the critical density): fl^ = 0.241 ±0.034. Ordi- 
nary matter parameter density (density of mass of ordi- 
nary atomic matter in the universe divided by the critical 
density): fit = 0.0416±0.001. Density parameter of dark 
energy (energy density of dark energy in the universe di- 
vided by the critical density): = 0.759 ± 0.034 ^. 

Even though of all successes of the ACDM, this model 
has several problems. Some of them are: The exotic 
weakly interacting particles proposed as dark matter par- 
ticles candidates are still undetected in the laboratory. 
The number of satellites in a galaxy such as the Milky 
Way is predicted to be an order of magnitude larger than 
is observed. Cuspy halo density proflles. The lack of 
evidence in the Milky Way for a major merger is hard 
to reconcile with the amount of accretion predicted by 
ACDM. With respect to the inclusion of the cosmologi- 
cal constant, the ratio of the vacuum energy density to 
the radiation energy density after inflation is 1 part in 
10^°°, a flne tuning coincidence. The cosmological con- 
stant has the wrong sing according to the string theorists 
who prefer a negative A instead of a positive A. The 
ACDM model predicts that large structures should form 
last and therefore should be young whereas observations 
tell us that the largest galaxies and clusters appear old. 
Several candidates have been proposed in the past that 
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pretend to substitute the role that the cosmological con- 
stant plays to accelerate the expansion of the universe. 
That the universe is expanding is supported mainly by 
the observations of the supernova project SNIa[^,!8!] . And 
we have lead to conclude that the universe is now dom- 
inated by an energy density with negative pressure and 
occupies about 70% of the universe in the present epoch. 
This energy is called generically the dark energy. Sev- 
eral models to explain this dark component has been 
proposed, and they may be classified accordingly to its 
equations of state as the following: quintessence dark 
energy 0, [13], phantom energy [ll|, [l^l, or the quintom 
cosmology paradigm (see also the review for even 
more details). Cosmological constant has as its equation 
of state one in which the pressure is the negative of the 
density. 

From the iV-body numerical simulations point of view 
we have two works, in which a detailed analysis of two 
scalar field possibilities has been explored. One is the 
coupled quintessence models [1^ and the other is the 
extended quintessence models [16|. The former considers 
a scalar field coupled minimally to the Ricci scalar and 
is coupled to the Lagrangian matter contribution. The 
later is of the type of a scalar-tensor theory in which the 
scalar field is introduced to model dark energy with a 
potential that is an inverse power law. The model we 
present in this work pretends to model the dark matter 
contribution in the large scale structure formation with 
a scalar field model that stems from the newtonian limit 
of a general scalar-tensor theory [TtI . 

Another problem that the ACDM problem is facing is 
the following. Almost a decade ago a bow shock in the 
merging cluster 1E0657-56, known as the Bullet Cluster, 
observed by satellite Chandra indicates that the subclus- 
ter -found by fisj- moving through this massive {lO^^h^^ 
Mq) main cluster creates a shock with a velocity as high 
as 4700 km s~^ [H, H^. A significant offset between the 
distribution of X-ray emission and the mass distribution 
has been observed [21, 22J, also indicating a high- velocity 
merger with gas stripped by ram pressure. Several au- 
thors have done detailed numerical non-cosmological sim- 
ulations (23l - [27| . One of the key input parameters for the 
simulations is to set the initial velocity of the subcluster, 
which is usually given at somewhere near the virial radius 
of the main cluster. [2^ have run cosmological N-body 
simulation using a large box (27 Gpc'^) to calculate 
the distribution of infall velocities of subclusters around 
massive main clusters. The infall velocity distribution 
were given at 1-3 i?200 -similar to the virial radius- and 
thus it gives the distribution of realistic initial velocities 
of subclusters just before collision. This distribution of 
infall velocities must be compared with the best initial 
velocity used by [l^] of 3000 km s~^ at about 2i?2oo to 
be in agreement with observations. [1^ have found that 
such a high infall velocity is incompatible with the pre- 
diction of the ACDM. 

Therefore, there are plenty of problems that the con- 
cordance model have to solve and finally, it does not tell 



us what is dark matter and dark energy. 

The program to study the large scale structure forma- 
tion should be to start with primordial initial conditions 
which means give the initial relevant fields, such as for 
example, density and velocity fields at the epoch of last 
scattering (~ z = 1100, the value of the redshift at that 
epoch, i.e., a photon emitted at that epoch is redshifted 
as 1 -|- z = Xr/Xe, with Ae the wavelength of the pho- 
ton when emitted, and is the wavelength of the same 
photon when observed; the expansion factor for a uni- 
verse with a flat geometry is related to the redshift as 
1 + z — 1/a). Then, evolve this initial condition using an 
A^-body scheme up to the present epoch (z = 0). 

Some questions we have to answer are: What arc the 
distribution of the LSS sizes? What is the amount of 
mass and its distribution at large scales? How are the 
voids distributed through the space? Are these voids 
devoid of any matter? How the LSS evolve with time? 
What is the DM equation of state? What is its role 
in the LSS formation processes and galactic dynamics? 
What are the implications of the observed LSS on the 
cosmological model of our universe? And on the struc- 
ture formation? And of course, what is the nature of the 
dark energy and matter? 

During the last decades there have been several pro- 
posals to explain DM, for example: Massive Compact 
Halo Objects (Machos), Weakly Interacting Massive Par- 
ticles (WIMPs) such as supersymmetric particle like the 
neutralino. Other models propose that there is no dark 
matter and use general relativity with an appropriate 
equation of state. Or we can use scalar fields, minimally 
or non minimally coupled to the geometry. 

In this work we are mainly concern with the problem 
of dark matter and its consequences on the large scale 
structure formation process. Our DM model is based on 
using a scalar field (SF) that is coupled non-minimally 
to the metric through the Ricci scalar in the Einstein 
field equations. A scalar field is the most simple field of 
nature. Nordstrom proposed a gravity theory by 1912, 
before Einstein [2^ . Scalar fields have been around for so 
many years since pioneering work of Jordan, Brans, and 
Dicke. Nowadays, they are considered as: (a) inflation 
mechanism; (b) the dark matter component of galaxies; 
(c) the quintessence field to explain dark energy; and so 
on. Therefore, is natural to consider dark matter mod- 
els based on modifications of Einstein's general relativity 
that include scalar fields. In this paper we will show 
results about the role this scalar field plays on the non- 
linear large scale structure formation of the universe. In 
particular, we will show how the power spectrums pre- 
dicted by this model compare with the power spectrums 
predicted by ACDM and the ones that come from obser- 
vations. 

So we organize our work in the following form: In the 
next section we present the general theory of a typical 
scalar-tensor theory (STT), i. e., a theory that general- 
izes Einstein's general relativity by including the contri- 
bution of a scalar field that couples non-minimally to the 
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Ricci scalar. In section Hill we show how the Friedmann 
equations becomes within a STT and present our model 
for the evolution of the universe expansion factor a{t). In 
section HV] we present the iV-body method which we will 
use to obtain the evolution of the large scale structures. 
Our results for an initial condition of the fields that is 
consistent with the observations are given in section [V] 
Finally, our conclusions are given in section IVTl 



II. GENERAL SCALAR-TENSOR THEORY 
AND ITS NEWTONIAN LIMIT 

The Lagrangian that gives us the Einstein equations 
of general relativity is 



for the metric g/^^i, and for the massive SF (p we have 



2oj 



[8TTT-uj'{d(l>f + <j)V' -2V] 



(6) 



where ()' = Here T^^ is the energy-momentum ten- 
sor with trace T, Ld{(t>) and V{(f>) are in general arbi- 
trary functions that govern kinetic and potential contri- 
bution of the SF. If in Lagrangian dU we set V{(f>) = 
we get the Bergmann- Wagoner theory. If we further set 
= constant the Jordan-Brans-Dicke theory is recov- 
ered. The gravitational constant is now contained in cj). 
Also, the potential contribution, V{(j)), provides mass to 
the SF, denoted here by msF- 



C = 



R 



(1) 



The Einstein field equations that are obtained from the 
above Lagrangian, in the limit of small velocities as com- 
pared with the speed of light and small forces, limit 
known as the Newtonian limit, give us the standard New- 
tonian potential due to a point particle of mass m 



N 



-G- 



(2) 



where G is the gravitational constant. What we intend 
to do in this work is to obtain the consequences in the 
LSS formation processes when we rise the constant G to 
a scalar field, 1/G — > (/). But we go beyond this approach 
and include in the Lagrangian two additional terms that 
depend on this field, a kinetic and potential terms. We 
will show, in particular, that the Newtonian limit of this 
theory gives for the Newtonian potential due to a mass 
m 



N 



1 



(3) 



i. e., the standard Newtonian potential is modified by an 
additional term that has the form of a Yukawa potential. 

Then, we start with the Lagrangian of a general scalar- 
tensor theory 



C 



IGtt 



+ Cm{9^.u) ■ (4) 



Here g^^i, is the metric, CMig^j.^) is the matter Lagrangian 
and ijj{(t>) and V^cj)) are arbitrary functions of the scalar 
field. The fact that we have a potential term V{(l>) tells 
us that we are dealing with a massive scalar field. Also, 
the first term in the brackets, (f>R, is the one that gives 
the name of non-minimally coupled scalar field. 

When we make the variations of the action, S = 
J d'^xC, with respect to the metric and the scalar field 
we obtain the Einstein field equations [1^ 



7? -In ff-i 



1 UJ 



, (5) 



A. Newtonian limit of a STT 

The study of large-scale structure formation in the 
universe is greatly simplified by the fact that a limit- 
ing approximation of general relativity, the Newtonian 
mechanics, applies in a region small compared to the 
Hubble length cH'^ {cH^^ ~ 3000/i~i Mpc, where c 
is the speed of light, = lOO/i km/s/Mpc, is Hubble's 
constant and h « (0.65 — 0.75)), and large compared to 
the Schwarzschild radii of any collapsed objects. The 
rest of the universe affect the region only through a tidal 
field. The length scale cH^^ is of the order of the largest 
scales currently accessible in cosmological observations 
and ^^(7^ ~ 10^°/i~^ yr characterizes the evolutionary 
time scale of the universe [3lj |. 

Therefore, in the present study, we need to consider 
the infiuence of SF in the limit of a static STT, and then 
we need to describe the theory in its Newtonian approxi- 
mation, that is, where gravity and the SF are weak (and 
time independent) and velocities of dark matter particles 
are non-relativistic. We expect to have small deviations 
of the SF around the background field, defined here as 
(0) and can be understood as the scalar field beyond all 
matter. Accordingly we assume that the SF oscillates 
around the constant background field 



and 



=(</))+ 



where t^^^ is the Minkowski metric. Then, Newtonian 
approximation gives [sO, l32l - [33 | 



Roo — TC^^hi^o ~ 



'■SFl 



G 



N 



4,71 p - 



1-t- a 



1 . 
2^^ 



(7) 
(8) 



we have set (0) — {l + a)/GN and a = 1/(3-1- 2cj). In the 
above expansion we have set the cosmological constant 
term equal to zero, since on small galactic scales its in- 
fiuence should be negligible. However, at cosmological 



4 



scales we do take into account the cosmological constant 
contribution, see below. 

Note that equation ([7]) _can be cast as a Poisson eqna- 
tionfor V= (l/2)(/ioo + 0/(0)), 



V'V' = 47r 



G 



N 



1 



(9) 



and the New Newtonian potential is given by $jv = 
(1/2) /loo ~ "0 ^ Above equation together 

with 



Sirap 



(10) 



form a Poisson-Helmholtz equation and gives 

1 Gat 7 



21 + a 

which represents the Newtonian limit of the STT with 
arbitrary potential V{(t)) and function a;(0) that where 
Taylor expanded around (0). The resulting equations 
are then distinguished by the constants Gat, a, and A — 
hp/mspc- Here ft,p is Planck's constant. 

The next step is to find solutions for this new New- 
tonian potential given a density profile, that is, to find 
the so-called potential-density pairs. General solutions 
to Eqs. © and pU)) can be found in terms of the cor- 
responding Green functions, and the new Newtonian po- 
tential is ISO, Sal 



G 



N 



1 + a 
Gn 



dr. 



1 



dv,^^ . +B.C. (11) 



The first term of Eq. pT|) . is the contribution of the 
usual Newtonian gravitation (without SF), while infor- 
mation about the SF is contained in the second term, 
that is, arising from the influence function determined 
by the modified Helmholtz Green function, where the 
coupling Lo (a) enters as part of a source factor. 

The potential of a single particle of mass m can be 
easily obtained from (|lip and is given by 



N 



Gn m 
1 + a r 



(l + ae-'^/^) (12) 



For local scales, r ^ A, deviations from the Newtonian 
theory are exponentially suppressed, and for r ^ A the 
Newtonian constant diminishes (augments) to Gn/O^ + 
a) for positive (negative) a. This means that equation 
fulfills all local tests of the Newtonian dynamics, and 
it is only constrained by experiments or tests on scales 
larger than -or of the order of- A, which in our case is of 
the order of galactic scales. In contrast, the potential in 
the form of equation ([3]) with the gravitational constant 
defined as usual does not fulfills the local tests of the 
Newtonian dynamics (36j . 

It is appropriate to give some additional details on 
the Newtonian limit for the Einstein equations with- 
out scalar fields (see [Slj). We are considering a small 



region compared to the Hubble length cH~^ but large 
compared to the Schwarzschild radii of any collapsed ob- 
ject. In this small region the metric tensor was written 
g^jy — rjf^i, -\- hf^i, where h^^, is small as compared to 
the Minkowski metric rjf^^. In this region Einstein's field 
equations are simple because the standard weak field lin- 
ear approximation applies. One finds. 



goo = + 2$ 



(13) 
(14) 



Then, the zero-zero component of the field equations for 
an ideal fluid with density p, pressure p, and velocity 
V <^ c becomes 



(15) 



For completeness the cosmological constant has been 
added. The geodesic equations, in the limit f <C c, 
h 1, are 



IF 



(16) 



Equations (ITSI) and are the standard equations of 
Newtonian mechanics, except that if there is an appre- 
ciable radiation background, one must take into account 
the active gravitational mass associated with the pres- 
sure, and of course if A ^ 0, there is the cosmic force 
Ar/3 between particles at separation r. 

Equations (ITSl) and (ITB)l apply to any observer out- 
side a singularity, though depending on the situation, 
the region within which these equations apply need not 
contain much matter. The region can be extended by 
giving the observer an acceleration gi to bring the ob- 
server to rest relative to distant matter, which adds the 
term g^r* to $, and then by patching together the results 
from neighboring observers. This works (the acceleration 
and potentials can be added) as long as relative velocities 
of observers and observed matter are ^ c and ^ 
(equation (jl4p ). For a region of size R containing mass 
M ~ pB? with density p roughly uniform, this second 
condition is 

GnpR^ < (17) 

In the Friedmann-Lemaitre models Hubble's constant is 

H^{GnpY'^ (18) 

If one assumes A is negligible and the density parameter 
^ 1, so equation ([T7| indicates 



R < cH-^ - 3000 Mpc - 10^** cm 



(19) 



That is, the region must be small compared to the Hubble 
length. Since the expansion velocity is u iJr, this 
condition also says v <^ c. 
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The Newtonian approximation can fail at much smaller 
R if the region includes a compact object like a neutron 
star or black hole, but one can deal with this by noting 
that at distances large compared to the Schwarzschild 
radius the object acts like an ordinary Newtonian point 
mass. It is speculated that in nuclei of galaxies there 
might be black holes as massive as 10^ Mq, Schwarzschild 
radius ~ 10^^ cm. If this is an upper limit, Newtonian 
mechanics is a good approximation over a substantial 
range of scales, lO^'* cm ^ r ^ 10^^ cm. 



B. Multipole expansion of the Poisson-Helmholtz 
equations 

The Poisson's Green function can be expanded in 
terms of the spherical harmonics, Yin{0, (p), 

' ' 1=0 n=~l > 

where r< is the smaller of |r| and jr^l, and r> is the 
larger of |r| and Ir^l and it allows us that the standard 
gravitational potential due to a distribution of mass p(j), 
without considering the boundary condition, can be writ- 
ten as 

^(r) = V-^'^ + ij'^"^ 

where (V'*-'^'') are the internal (external) multipole 
expansion of -0, 

00 / 



and its force is 







1=0 n=-l 
00 / 

EE 

1=0 n=-l 



47r (e)Yin{e,if) 



21 + f^'" r'+i 



Here, the coefficients of the expansions ■0^'-' and 
known as internal and external multipoles, respectively, 
are given by 



47r / dr'^y,;(0',^')p(r'): 

V(r<r') 

47r/ dr%l{e\v')r" p{r') . 

V(r>r') 



The integrals are done in a region V where r < r' for 
the internal multipoles and in a region V where r > 
for the external multipoles. They have the property 

(20) 



We may write expansions above in cartesian coordi- 
nates up to quadrupoles. For the internal multipole ex- 
pansion we have 



r p 



1 



-r Q 



(0 



(21) 



F«=p«-l-Q«-r, 



where 



Jv(r<r') 



Jv(T<r') r'-^ 



Jv(r<r') ^ 

For the external multipoles we have 

_ M(^) r • pf'^) 1 r • Q^'^) • r 



and its force is 



r(e) _ 



j^3 ^3 



where 



]V&^ = f dv'piv') , 

Jv{r>r') 



pf^^ / dr'x'Av')., 

Jv{r>r') 



V{r>r') 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



dr' {3x',x'j-r'^S,j)p{r'). (30) 



The external multipoles have the usual meaning, i.e., 
M^*^^ is the mass, p^^^ is the dipole moment, and Q'-'^' is 
the traceless quadrupole tensor, of the volume V{r > r'). 
We may atach to the internal multipoles similar mean- 
ing, i.e., M'^*^ is the internal "mass", p*^*^ is the inter- 
nal "dipole" moment, and Q^*-* is the traceless internal 
"quadrupole" tensor, of the volume V(r < r'). 

In the case of the scalar field, with the expansion 

exp(-TO|r - r^l) • / m / n 
— i = 47rTO } ^ y ^ iiirnr^jkiimry) 



the contribution of the scalar field to the Newtonian grav- 
itational potential can be written as 



(r) 



6 



where, for simplicity of notation, we are using m = 

msF = hp/{c\) and 



OO I 



1 _ ^ <- 



1=0 n=-l 



ii{x) and ki{x) are the modified spherical Bessel func- 
tions. 

We have defined the multipoles for the scalar field as 
€i=^ I dr' ^'-y^ {mr'y+'ki{mr') p(r'), 

Jv{r<r') 



9' il^^l^^n 

IV{r>r') ' {"^rj 

They, also, have the property 



= (-i)"(C)* 



(31) 



The above expansions of SF contribution to the New- 
tonian potential can be written in cartesian coordinates. 
The internal multipole expansion of the SF contribution, 
up to quadrupoles is 



2a mr 



2 {mry 



(32) 



and its force is 



= -m^!l(^MWr - 3^p« (33) 



mr mr 
-3m^-p (p^ • r)r - 5-- r^Q^ ' ■ r 



+5^[5^-^](r-Q«-r)r, 
2r^ (mr)^ mr 



where 



mW = /" dr'(mr')A;o(mr')^p(r'), (34) 
#=/ dr'(mr')'fci(mr')a;;4Mr'), (35) 



(i) 



V(r<r') 



dr' (mr')^fc2(mr')(3a;^a;^ -r'2<5ij)^p(r') . 

(36) 

In the exterior region the SF multipole contribution to 

the potential is 



2a 



=mrko{mr) h 3(mr)^A;i(mr)- 

1 , ,3, , ,r •Q(^) -r 
+o-(mr) K2[mr) g 



and its force is 

— F^^ = (mr) ki{mr) — §— r— 3(mr) fci(mr) — ^ 



2a 



-|-3 [(mr)fco(mr) 

-|-(2 -f mr){mr)'^ ki{mr)] 
Q(«=) • r 



(p(e) . r) 



(38) 



where 



— 5(mr)^fc2(mr)- g 
5 

+ 2 [3(™^)fco(™'') + 3(rrtr)^A;i(m,r) 

(r • Q(^) • r) 
-I- (3-hmr)(mr)^fc2(mr)] — ^r. 



M(«) = / dr'io{mr')p{r') , (39) 

Jv{r>r') 



Jv{r>r' 



, ilimr') , 



dr' 

/) mr 



7^^iP(r'), 



(40) 



i2(mr') , , 



'5? - / '^r' j^i^x'.x', - r'%Mr') ■ (41) 

Jv{r>r') ) 

In the limit when m — >■ we recover the standard New- 
tonian potential and force expressions. 

Up to here the formulation is general, i.e., mass distri- 
bution may have any symmetry or none at all. In order to 
take advantage of the symmetry of the spherical harmon- 
ics, the mass distribution must be spherically symmetric. 



III. COSMOLOGICAL EVOLUTION 
EQUATIONS USING A STATIC STT 

To simulate cosmological systems, the expansion of the 
universe has to be taken into account. Also, to determine 
the nature of the cosmological model we need to deter- 
mine the composition of the universe, i. e., we need to 
give the values of ili = Pi/pc, with = SH"^ /SttGn, 
for each component z, taking into account in this way all 
forms of energy densities that exist at present. If a par- 
ticular kind of energy density is described by an equation 
of state of the form p = wp, where p is the pressure and w 
is a constant, then the equation for energy conservation 
in an expanding background, d{pa^) = —pd{a^), can be 
integrated to give p cx a^^'-^^'^\ 

Then, the Friedmann equation for the expansion factor 
a{t) is written as 



!i=^oE".( 



(42) 



(37) 

where Wi characterizes equation of state of species i. 
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The most familiar forms of energy densities are those 
due to pressureless matter with Wi = (that is, nonrel- 
ativistic matter with rest-mass-energy density pc^ domi- 
nating over the kinetic-energy density pv'^ /2) and radia- 
tion with Wi = 1/3. The density parameter contributed 
today by visible, nonrelativistic, baryonic matter in the 
universe is fis ~ (0.04 — 0.05) and the density parameter 
that is due to radiation is fin 2 x 10~^. 

In this work we will consider a model with only two 
energy density contribution. One which is a pressureless 
and nonbaryonic dark matter with Udm ~ 0.3 that does 
not couple with radiation. Although in the numerical 
simulations we may include in ^Idm the baryonic matter. 
Other, that will be a cosmological constant contribution 
r^A ~ 0.7 with and equation of state p = —p. The above 
equation for a{t) becomes 



\ a 



(43) 



The above discussion gives us the standard cosmo- 
logical model with cosmological constant, i. c., ACDM 
model. 

In the framework of a scalar-tensor theory the cosmol- 
ogy is given as follows. If we use the Friedmann metric 



dr^ 



1 — Kr^ 



r'^{de 



2 ' sin2 



(44) 

in the time-time component of the Einstein field equa- 
tions (Hamiltonian constraint) gives 



8^ ^0 ^(0) ( 
3</>^ 4> & \4> 



K 



V 

60 



(45) 



while the equation for the scalar field is 



3H 



2uj + 3 



1 



dV 



2uj + 3 

and the equation of the fluid is 



Stt (p - 3P) - 0— + 2V 



p + 3H{p + P)=0 



(46) 



(47) 



The cosmological evolution of the initial perturbed 
fields should be computed using the above equations for 
the expansion factor. However, here, we will employ a 
cosmological model with a static SF which is consistent 
with the Newtonian limit given above. Thus, the scale 
factor, a(t), is given by the following Friedman model. 



n 



DMO 



1 



1 



DMO 



AO 



1 + a 

' (48) 

where H = a/a, floMo and Q\o Siie the matter and en- 
ergy density evaluated at present, respectively. The de- 
nominator (1 -I- a) appears in Eq. (|48p due to the fact 



CMB Anisotropies 




Multipole Moment / 



FIG. 1. Cosmic Microwave Background anisotropies. With 
small dash symbol are the WMAP 7 years experimental ob- 
servations [43]. We have plotted several cases in which the 
Newton constant has been modified. The standard case is 
shown with a small dash line. With a long dash Une is shown 
an efi^ective gravitational constant, Ge// = 0.5G]v. With a 
dotted Une is shown the case with Ge// — 1.5Gn, and with a 
dash-dotted line is shown the case with Ge// ~ 2Gn- 



that we have defined fli = p/pc and pc — ^H'^/SttGn. 
Then, whenever appears the gravitational constant Gn, 
we replaces it by Gat/ (1 + a). 

We notice that the source of the cosmic evolution is de- 
viated by the term 1 -I- a when compared to the standard 
Friedmann-Lemaitre model. Therefore, it is convenient 



to define a new density parameter by fl. 



(a) 



n,/{i + a). 



This new density parameter is such that i^^j// + ^a''' = 
1, which implies a flat universe, and this shall be as- 
sumed in the following computations, where we consider 

in^^\n^j^^) = (0.3,0.7). 

For positive values of a, a flat cosmological model de- 
mands to have a factor (1 -I- a) more energy content {Qm 
and I^a) than in standard cosmology. On the other hand, 
for negative values of a one needs a factor (l + a) less 
and Q\ to have a flat universe. To be consistent with the 
CMB spectrum and structure formation numerical ex- 
periments, cosmological constraints must be applied on 
a in order for it to be within the range (—1,1) [38l - l4ll |. In 
figured] we show the effect of different values of the grav- 
itational constant Gat on the anisotripies of the CMB. 



IV. VLASOV-POISSON-HELMHOLTZ 
EQUATIONS AND THE TV-BODY METHOD 

The Vlasov-Poisson equation in an expanding uni- 
verse describes the evolution of the six-dimensional, one- 
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particle distribution function, /(x, p). The Vlasov equa- 
tion is given by, 



at ma'^ ox op 







(49) 



where x is the comoving coordinate, p = ma^x., m is the 
particle mass, and $Ar is the self-consistent gravitational 
potential given by the Poisson equation. 



WHn{x) = AttGn a" [p(x) - pb{t)] 



(50) 



where is the background mass density. The Vlasov- 
Poisson system, Eqs. ((55)) and (1501 form the Vlasov- 
Poisson equation, constitutes a collisionless, mean-field 
approximation to the evolution of the full iV-body distri- 
bution. 

An TV-body code attempts to solve the Vlasov-Poisson 
system of equations by representing the one-particle dis- 
tribution function as 



N 



/(x, p) = ^ ^(x - Xi)5{p - Pi) 



(51) 



Substitution of (|5T|) in the Vlasov-Poisson system of 
equations yields the exact Newton's equations for a sys- 
tem of N gravitating particles (See ^3 details). 



1 X ^ mj{xi 
2Hxi — tGjv > — i 



(52) 



where the sum includes all periodic images of parti cle j 
and numerically is done using the Ewald method, see [4J] . 
It is important to keep in mind, however, that we are not 
really interested in solving the exact TV-body problem for 
a finite number of particles, N. The actual problem of 
interest is the exact A^-body problem in the fluid limit, i. 
e., as A'^ — >■ oo. For this reason, one important aspect of 
the numerical fidelity of A^-body codes lies in controlling 
errors due to the discrete nature of the particle represen- 
tation of /(r, p). 

In the Newtonian limit of STT of gravity to describe 
the evolution of the six-dimensional, one-particle distri- 
bution function, /(x, p,t) we need to solve the Vlasov- 
Poisson-Helmholtz equation in an expanding universe 
[45l - l49| . The Vlasov equation is given by 



5/ 
dt 



^ — -mV$^^(x)•|^ = 



dp 



(53) 



where = (l/2)ft,oo 
Poisson equation 

VV(x) 



tp — ^j^(t> with ip satisfying 
47rGNa^[pi^)~Pbit)] (54) 



and (j) satisfying the Helmholtz equation 

\7^4>-X-^4> = -8Traa^[p{x)-pb{t)] , (55) 

Above equations form the Vlasov-Poisson- Helmholtz sys- 
tem of equations, constitutes a collisionless, mean-field 



approximation to the evolution of the full A^-body dis- 
tribution in the framework of the Newtonian limit of a 
scalar-tensor theory. 

Using the representation of the one-particle distribu- 
tion function (j5ip the Newtonian motion equation for a 
particle i, is written as [46| 

^.+2hx.^-i,^j:^¥^ 

a-^ I + a ^ \xi — X, P 
xi^si?(|x., - Xj|,Q!, A) (56) 

where the sum includes all periodic images of particle j, 
and FspiT, a. A) is 



SF 



(r,a,A) = 1 + a (l + J) e^"^/^ 



(57) 



which, for small distances compared to A, is Fsp(r < 
A, a. A) w 1 + a [l + j) and, for long distances, is 
Fspif > A, a. A) ~ 1, as in Newtonian physics. 

The function Fsrir, a, A) is shown in Fig. [2]for several 
values of parameter a. The horizontal line at Fsf = 1 
is for the standard Newtonian case {a = 0). The long 
dash line above the horizontal standard line is for a = 1, 
the medium dash line is for a = 1/2 and the tiny dash 
line is for a = 1/4. For negative values of a the lines are 
below the standard horizontal case. Long dash line is for 
a = —1/2 and medium dash line is for a — —1/4. We 
should notice that even though A = 5 Mpc/h, Fsf gives 
an important contribution to the force between particles 
for r > A. 



FIG. 2. Function FsF{r,a,\) for A = 5 Mpc/h and several 
values of parameter a. a = 0, i. e., the standard Newtonian 
case is the horizontal line. The long dash line above the hori- 
zontal standard line is for a — 1. The medium dash line is for 
a — 1/2 and the tiny dash line is for a — 1/4. For negative 
values of a the lines are below the standard horizontal case. 
Long dash line is for a = —1/2 and the medium dash line is 
for a = -1/4. 



V. RESULTS 

In this section, we present results of cosmological sim- 
ulations of a ACDM universe with and without SF con- 
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tribution in order to study the large scale structure for- 
mation. 

We have used the standard Zel'dovich approximation 
^ to provide the initial 256"^ particles displacement off 
a uniform grid and to assign their initial velocities in a 
256 Mpc box[50|- In this approximation the comov- 
ing and the lagrangian coordinates are related by 



«2y"fo|k|(t)S|k|(q) 



k 



,2 \ ^ 



k 



&ik|W 



&ik|W 



(58) 



S|k|(q), (59) 



where the displacement vector S is related to the veloc- 
ity potential $ and the power spectrum of fluctuations 
P(|fc|): 



S|k|(q) = V,$|k|(q), 



(60) 



*|k| = X! ^"'^ • q) + &k sin(k • q)) (61) 



o 



(xlO^) 



(a) Without SF 




0.5 1.0 1.5 2.0 2.5, 
X (kpc/h) (xlO ) 



where a and b are gaussian random numbers with the 
mean zero and dispersion tr^ — P{k)/k'^, 



Ok 



G(0,1), 



&k 



G(0,1) (62) 



where G(0, 1) is a gaussian number with mean zero and 
dispersion 1. 

The parameter az together with the power spectrum 
P{k), define the normalization of the fluctuations. The 
initial power spectrum was generated using the fitting 
formula by ^51,] for the transfer function. This formula is 
a slight variation of the common BBKS fit . 

Therefore, for the standard ACDM we have for the 
initial condition that the starting redshift is = 50 
and we choose the following cosmology: ftoM = 0.314 
(where ^dm includes cold dark matter and baryons), 
flB ^ 0.044, = 0.686, Hq = 71 km/s/Mpc, 

as = 0.84, and n = 0.99. Particle masses are in the order 
of 1.0 X 10^" Mq. The individual softening length was 
50 kpc/h. This choice of softening length is consistent 
with the mass resolution set by the number of particles. 
All these values are in concordance with measurements 
of cosmological parameters derived from the seven-year 
data of the WMAP |6|. The initial condition -called the 
big box case- is in the Cosmic Data Bank web page: 



FIG. 3. (a) LSS formation process without the SF contribu- 
tion. 



field a flat universe is obtained if D, 



1, with 



(http : //tSweb . lanl .gov/people/heitinciiin/test3 .html) 

See [50[] for more details. 

Because the visible component is the smaller one and 
given our interest to test the consequences of including 
a SF contribution to the evolution equations, our model 
excludes gas particles, but all its mass has been added to 
the dark matter. We restrict the values of a to the inter- 
val (-1,1) (see [38-41]) and use A = 1,5,10,20 Mpc/h, 
since these values sweep the scale lengths present in the 
simulations. 

In all the simulations we have done, with or without 
the scalar field contribution, we demand that the cos- 
mological model be fiat. In this model with a scalar 



n\n' = ^„i/{l + a) and fl'-j^' = n^/{l + a). Then, for 
positive values of a we need a factor (l + a) more energy 
content {Qm and Q\) that in the standard cosmology. 
Whereas for negative values of a we need a factor of 
(1-1- a) less flm and 17a to have a flat universe. With this 
recipe and for a given value of a we modified the above 
ACDM initial condition accordingly. 

In Fig. [3] we show how the above initial conditions 
evolve and give us the LSS formation process without SF 
and with SF. In (a) without SF. In (b) with SF: a = 1 
and A = 1 Mpc/h. In (c) with SF: a = 1 and A = 5 
Mpc/h. In (d) with SF: a = 1 and A = 10 Mpc/h. In 
(e) with SF: a 1 and A = 20 Mpc/h. In (f) With SF: 
a = -1/2 and A 5 Mpc/h. In (g) With SF: a = -1/4 
and A = 5 Mpc/h. 

To study the structure formation in the universe we 
follow the evolution of the overdensity 

,5(x) = 1 

Pa 

where po is the average density over a volume V and x 
is the comoving distance related to the physical density 
by r = a(t)x. In the linear regime S <^ 1. 

The correlation function tell us how 6 is correlated in 
two nearby points x' and x' + x, 

e(x)^(J(x')<5(x'+x)) 

and the power spectrum is the Fourier transform of the 
correlation function. 
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FIG. 3. (b)-(c) LSS formation process with the SF contri- 
bution, (b) Q = 1 and A = 1 Mpc/ft. (c) a = 1 and A = 5 
Mpc/h. 



FIG. 3. (d)-(e) LSS formation process with the SF contribu- 
tion, (d) Q = 1 and A = 10 Mpc/h. (e) a = 1 and A = 20 
Mpc/h. 



In Fig. m (a) we show the evolution of the power 
spectrum for the big box L — 256 Mpc/h without SF, 
for several vahies of the redshift z. We have used the 
POWMES code to compute the matter power spectrum 
[53| . In this figure we can appreciate how it is forming 
structures, upper curve of the power spectrum for the 
present epoch, z = 0. Fig. |3] (b) shows the power spec- 
trum for the same values as in Fig. [31 Continuous hue is 
without SF. Upper curve (dashed line) is with SF: a = 1 
and A = 5 Mpc/h. Lower curve (dashed-dotted line) is 
with SF: a = -1/2 and A = 5 Mpc/h. The curve that 
is just below the continuous line (dotted line) is with 
SF: a = —1/4 and A = 5 Mpc/h. More greater values 



of the power spectrum means more structure formation. 
Therefore, the inclusion of a SF modifies the structure 
formation process. Depending of the values of its pa- 
rameters, a and A we can obtain more structure or less 
structure at the present epoch, z = 0. 

In our results as shown in Figs. [3] (c), (f) and (g), 
and m we have used a fix value of A = 5 Mpc/h. For 
a given A the role of SF parameter a on the structure 
formation can be inferred by looking at equations (j56p 
and ((57)) and Fig. [2] where we show the behavior of Fsf 
as a function of distance for several values of a. The 
factor Fsf augments (diminishes) for positive (negative) 
values of a for small distances compared to A, resulting 
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FIG. 3. (f)-(g) LSS formation process with the SF contribu- 
tion, (f) a = -1/2 and A = 5 Mpc//i. (g) a = -1/4 and 
A = 5 Mpc/h. 



in more (less) structure formation for positive (negative) 
values of a compared to the ACDM model. In the case 
of the upper curve in Fig. |4] (b), for r ^ A, the effective 
gravitational pull has been augmented by a factor of 2, in 
contrast to the case shown with the lower line in Fig. |4] 
(b) where it has been diminished by a factor of 1/2. That 
is why we observe for r < A more structure formation in 
the case of the upper curve in Fig. |l](b) and lesser in the 
case of the lower curve in the same figure. The effect is 
then, for a increasing positive a, to speed up the growth 
of perturbations, then the halos, and then the clusters, 
whereas negative values of a (a -H- —1) tend to slow 
down the growth. We also observe that for the large scale 
regimen of our simulations (fc < 0.1 /i/Mpc) they tend 



FIG. 4. (a) Evolution of the power spectrum for the case 
which there is no SF. (b) The power spectrum including SF 
for several values of parameter a and with A = 5 Mpc/ft. 
Experimental data are from galaxies in the Sloan Digital Sky 
Survey (SPSS) catalog and from an analysis of the Lyman-a 
forest ^41 



to predict almost the same structure formation. From 
comparison with the experimental results, we see that the 
ACDM agrees well with SDSS observations, but predicts 
more structure formation than observations show in the 
Lyman-a forest power spectrum. In general the more 
favored model is the model with SF with a = —1/2 and 
A = 5 Mpc/h. 
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VI. CONCLUSIONS AND FINAL COMMENTS 

We have used a general, static STT, that is compatible 
with local observations by the appropriate definition of 
the background field constant, i. e., ((/)) = [\ + a)/GN^ to 
study the LSS formation process. The initial condition 
for the several cases (different values of parameter a) 
was built in such a way that the geometry of the model 
universe were flat. Quantitatively, this demands that our 
models have f2/(l + a) = 1 and this changes the amount 
of dark matter and energy of the models in order to have 
a flat cosmology. 

Using the resulting modified dynamical equations, we 
have studied the LSS formation process of a ACDM uni- 
verse. We varied the amplitude and sign of the strength 
of the SF (parameter a) in the interval (—1, 1) and per- 
formed several SD-simulations with the same initial con- 
ditions. From our simulations we have found that the 
inclusion of the SF changes the local dynamical prop- 
erties of the most massive groups, however, the overall 
structure is very similar, as can be seen in Fig. [31 

The general gravitational effect is that the interaction 
between dark matter particles given by the potential 
(see equation changes by a factor Fsfi equation 

(|57p . in comparison with the purely Newtonian case. 
Thus, for a > the growth of structures speeds up in 
comparison with the Newtonian case (without SF). For 
the a < case the effect is to diminish the formation of 
structures. 

It is important to note that particles in our model are 
gravitating particles and that the SF acts as a mecha- 
nism that modifies gravity. The effective mass of the SF 
{msF — 1/A) only sets an interaction length scale for the 
Yukawa term. 

In this work we only varied the amplitude of the SF 
-parameter a- leaving the scale length. A, of the SF un- 



changed. However, in other studies we have done [48| 
we have found that increasing A enhances the structure 
formation process for a positive, and decreasing A makes 
the structures grow at a slower rate. 

We have computed the mass power spectrum in or- 
der to study the LSS formation process. The theoretical 
scheme we have used is compatible with local observa- 
tions because we have defined the background field con- 
stant < (f) >— + a) or equivalently that the lo- 
cal gravitational constant is given by (1 -|- a){4>)~^, in- 
stead of being given by {(f>)^^- A direct consequence 
of the approach is that the amount of matter (energy) 
has to be increased for positive values of a and dimin- 
ished for negative values of a with respect to the stan- 
dard ACDM model in order to have a flat cosmologi- 
cal model. Quantitatively, our model demands to have 
ri/(l + a) = 1 and this changes the amount of dark 
matter and energy of the model for a flat cosmological 
model, as assumed. The general gravitational effect is 
that the interaction including the SF changes by a factor 
FspifyCt, X) « 1 + a (l + for r < A in comparison 
with the Newtonian case. Thus, for a > the growth 
of structures speeds up in comparison with the Newto- 
nian case. For the a < case the effect is to diminish 
the formation of structures. For r > A the dynamics is 
essentially Newtonian. 

Comparison with the power spectrums from galaxies in 
the SDSS catalog and that inferred from Lyman-a for- 
est observations tell us that ACDM predicts more struc- 
ture formation in the regime of smaller scales {k > 0.4 
/i/Mpc). Whereas, the model with SF with a = —1/2 
and A = 5 Mpc/h, follows the general trend of the ob- 
served power spectrum. In this way we are able to con- 
struct a model that predicts the observed structure for- 
mation in the regime of small scales, with lesser number 
of halo satellites than the ACDM model. 
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